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ABSTRACT 

We conducted a very long baseline interferometric (VLBI) observation of 5 
nearby pulsars and did not resolve the scattering disks of any of these sources. 
Using our upper limits on the angular diameters of these scattering disks 
and published values of the broadening times and proper motion velocities, 
we constrain the possible distributions of scattering material. The material 
responsible for scattering these sources is neither uniformly distributed nor 
concentrated at the surface of the Local Bubble. We argue that these pulsars 
themselves influence their environments to produce this scattering material. 

Subject headings: turbulence, scattering - ISM: structure - pulsars: general - 
techniques: interferometric 

1. INTRODUCTION 

Fluctuations in the interstellar electron density scatter radio waves and cause angular 
broadening of radio sources. The angular sizes of these broadened sources contain 
information about the material responsible for scattering. Rao & Ananthakrishnan (1984) 
found that the angular diameters of low-latitude (|6| < 10°) radio sources tend to be larger 
in the direction of the galactic center. Dennison et al (1984) found large variations in the 
angular diameters of low-latitude extragalactic radio sources along nearby lines of sight. 
Some low-latitude sources are scattered by HII regions (Moran et al. 1990, Molnar et 
al. 1995) or supernova remnants (Spangler et al. 1986). These results suggest that the 
scattering material in the disk of our galaxy contains both an extended component and 
localized features. 

Pulsars are weak radio sources, but permit observation of the temporal broadening of 
their pulses as well as angular broadening. Cordes, Weisberg, & Boriakoff (1985) found 
strong variations in the temporal broadening of pulsars along nearby lines of sight. Gwinn, 
Bartel, & Cordes (1993) have compared the angular broadening and temporal broadening of 
pulsars to distinguish between different distributions of scattering material. They conclude 
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that the scattering properties of several distant pulsars are consistent with a uniform 
distribution of scattering material. The Vela pulsar is a notable exception; this young 
pulsar is scattered by its supernova remnant (Desai et al. 1992). 

Because of the small scale height of the galactic disk, radio sources at high galactic 
latitudes > 10°) are scattered by the local interstellar medium (ISM). Readhead 
& Hewish (1972) measured the angular diameters of high-latitude extragalactic radio 
sources and found that their angular sizes decreased smoothly with galactic latitude. They 
suggested that the scattering material is uniformly distributed above the galactic plane with 
a scale height of ~ 500 pc. In contrast, Hajivassiliou (1993) found in a similar survey that 
a model in which all scattering occurs at the edge of the Local Bubble could explain the 
measured distribution of angular diameters at high latitude. The Local Bubble is a cavity 
of hot coronal gas roughly 30-200 pc in extent that surrounds the earth (Cox & Reynolds 
1987). Nearby pulsars tend to have stronger refractive scattering on shorter timescales 
(Gupta, Rickett, & Coles, 1994) and more instances of multiple imaging (Gupta, Rickett, & 
Lyne 1994) than more distant pulsars. This may indicate that the distribution of material 
responsible for scattering differs in the two cases. One possibility is that nearby pulsars are 
scattered by the material at the surface of the Local Bubble, while more distant pulsars are 
scattered by a uniform medium. In this paper we describe a measurement to determine the 
location of material responsible for scattering nearby pulsars. 



2. INTERFEROMETRIC MEASUREMENTS OF ANGULAR 

BROADENING 

The large intensity variations of pulsars in both time and frequency result from 
scattering by electron density fluctuations in the ISM (Scheuer 1968). Fluctuations in the 
interstellar electron density produce fluctuaions in the index of refraction, thereby inducing 
variations in the relative phases of light rays propagating along different paths. The electric 
field measured at an antenna will vary as the phasor sum of these rays randomly coheres 
and decoheres. Rays that cohere at frequency v and time t remain coherent over a range of 
frequencies Az/ about v and over a range of times t; ss about t. Here Az/ is the scintillation 
bandwidth and ti SS is the scintillation timescale. An intensity maximum of bandwidth Av 
and duration t iss is called a "scintle" . 

An interferometer measures the visibility, which is the complex product of the electric 
fields at two antennas. Because the electric field at each antenna fluctuates, the visibility 
V(x, x + B;u,t) will be a function of the location of the two antennas, the frequency uj, 
and the time t. Here x and x + B are the positions of the antennas projected into the plane 
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perpendicular to the line of sight to the source. Averaging this visibility over position x 
gives 
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(V(x, x + B; to, t)) = V G exp 
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were V D is the average intensity measured at a single antenna and <p(x) is the phase of the 
electric field at position x. The equality follows from Mercier's theorem. We expect visibility 
fluctuations to be a stationary random process, so that averaging over x is equivalent to 
averaging over many scintles in time or frequency. The angular brackets ( ) indicate both 
kinds of averaging. We denote this averaged visibility as V{B). 

The average visibility V(B) is also the two-dimensional Fourier transform of the sky 
brightness distribution I(q), where B and q are conjugate variables. Expressions for V(B) 
have been derived for several different power spectra of electron density fluctuations; we use 
the results of Coles et al. (1987). For a power- law spectrum of density fluctuations with 
index a and inner scale L 



V(B) = V Q exp 
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where $h characterizes the angular diameter of the average scattering disk and A = 92 cm is 
our observing wavelength. The index a vis = a if \B\ > k, and « vis = 4 if \B\ < U. Estimates 
for the inner scale range from 10 2 km (Spangler & Gwinn 1990, Molnar et al. 1995) to 10 6 
km (Coles et al. 1987). We take a vis = 4 in this paper, corresponding to a Gaussian sky 
brightness distribution. In this case 9 H is the full width at half maximum (FWHM) angular 
diameter of the average scattering disk. Comparison of Equations 1 and 2 yields Oh'- 
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where D^B) = ((4>(x) — 4>{x + B)) 2 ) is the phase structure function. 



3. OBSERVABLE CONSEQUENCES OF SCATTERING 

A comparison of angular broadening and temporal broadening provides information on 
the distribution of scattering material. Expressions for the rms scattering angle 9 and the 
mean broadening time t of rays propagating through a distribution of scattering material 
have been derived by Alcock & Hatchett (1978) and Blandford & Narayan (1985). 
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Here z is a coordinate from the source at z = to the observer at z = D, ip(z) is the 
scattering rate per unit length, and c is the speed of light. The rms angular diameter of 
the average scattering disk is 9, and is related to the FWHM angular diameter 9h by 



® h = y 4 In 2 $ if a = 4. The mean broadening time t differs from the temporal broadening 
r by a factor close to unity (Lee & Jokipii 1975), which we neglect. If we assume a model 
of the scattering material ip(z) with one or two parameters, we can use measurements of 9h 
and r to derive the values of these parameters. 

In this paper we consider two special cases. If the scattering material is uniformly 
distributed along the line of sight, ip{z) is a constant and the FWHM angular diameter 

1 /2 

of the average scattering disk is 9 U = (161n2 ar/D) . If the scattering material is 
concentrated in a thin screen located a distance d s from earth, ip(z) oc 5(z — D + d s ) and the 
FWHM angular diameter of the average scattering disk is 9 S = (81n2cr (D — d s ) / Dd s ) 1 ^ 2 . 
We expect 9 H = 9 U for a pulsar scattered by a uniform medium. If the pulsar is scattered 
by a thin screen, we can use 9h and r to compute the distance d s to the screen. 

If the scattering material is concentrated in a thin screen, measurement of the proper 
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motion velocity V pm of a pulsar provides an estimate of the velocity V SCT of the scattering 
screen. Gupta et al. (1994) have discussed this relationship, which in our notation becomes 



Here all velocities are measured relative to the sun and we have assumed that the earth's 
orbital velocity is negligable. Gupta (1995) showed that the scattering properties and 
proper motion velocities of most pulsars are consistent with a screen halfway to the pulsar 
if V SCT = 0. By measuring the distance d s to the screen, we can estimate the velocity V^ cr of 
the scattering material. 



We attempted to resolve the scattering disks of 16 nearby pulsars with an 
intercontinental VLBI array. Our array consisted of telescopes at Jodrell Bank (76 m), 
Westerbork (phased array), Noto (30 m), and four 25 meter telescopes from the US Very 
Long Baseline Array of the National Radio Astronomy Observatory^]: Hancock, Owens 



1 NRAO is operated by the Associated Univesities Inc., under contract with the National Science 
Foundation. 
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OBSERVATIONS AND DATA REDUCTION 
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Valley, Brewster, and St. Croix. We observerd for 12 hours on 1994 May 25-26, recording 
left-circular polarization from 323 to 337 MHz with the Mark IIIA VLBI recording system. 
The observations were broken into 13 minute scans, and alternated between the pulsars 
and a set of extragalactic radio sources used for calibration. We detected only 2 calibration 
sources, 1345+125 and 2230+114, on intercontinental baselines. 

The tapes were correlated with the Mark IIIA correlator at the Max-Planck Institut 
fiir Radioastronomie in Bonn, Germany. We integrated the data over 2.5 second intervals 
and correlated over 24 lags, yielding 0.167 MHz frequency resolution. This is sufficient to 
resolve typical scintillation bandwidths of nearby pulsars at 326 MHz (see Table 1). Pulsar 
gating allowed us to correlate only during the pulse, yielding an improvement in the signal 
to noise ratio by a factor of 2-5 (Gwinn et al. 1986). The data from Hancock were discarded 
due to poor recording quality. 

We used Haystack Observatory's fringe fitting program FOURFIT to recover the 
cross-correlation functions from the Mark IIIA correlator output files. These were Fourier 
transformed to give the visibility V(x, x + B;u>,t). 

The visibility data were degraded by the fractional bitshift effect. This effect arises 
from the inability of a lag correlator to adjust the digitized data streams continuously 
to compensate for the variation of the geometrical delay r g (t). The sawtooth error in 
delay Ar g (t) causes a phase error 0fb s (c<j,t) = LuAr g (t). This error nearly averages out for 
unmodulated sources, but we found that beating between 0fb s (w,t) and the pulse gate 
introduced errors in the phases of as much as 0.5 radians. We computed the corrections 
to the visibility phases by reconstructing the sawtooth delay error from the Mark IIIA 
correlator output and matching it to the pulsar ephemeris. 

For each baseline in each scan we removed an instrumental model from the visibility 
phase. This model corrected for the phase errors introduced by the differing electronic 
pathlengths of the seven 2 MHz baseband converters, the geometrical delay, and the 
passband curvature within the baseband converters. Our phase correction for the 
geometrical delay contained 3 contributions: the multiband delay, the fringe rate, and the 
multiband acceleration. This last correction is not normally removed in interferometry, 
but was found to introduce errors of up to 0.5 radians on some baselines. The phase 
passband curvature of the baseband converters, pb (u;), was computed for each baseline by 
time averaging the visibility phases over the entire data set. These corrections were < 0.2 
radians. The resulting visibility phases were stable to 0.1 radians for our brightest sources. 
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Fig. 1. — The visibility of the pulsar B1919+21 plotted for the 7812 km Owens Valley - 
Westerbork baseline. On the left the visibility amplitude is plotted vs. time for 5 adjacent 
frequency channels out of 85 observed. The amplitude shows variations in both time and 
frequency caused by scattering in the ISM. The visibility phase is plotted vs. time on the 
right, and shows no such variations. This indicates that the scattering disk of the pulsar is 
unresolved on this baseline. The small fluctuations in the phase are consistent with noise. 



5. RESULTS 



Although we detected fringes on all but one of the 16 pulsars we observed, only the 5 
listed in Table 1 were bright enough to resolve individual scintles. All pulsars and quasars 
showed strong fluctuations of visibility phase with a bandwidth greater than 14 MHz. 
Scattering by the ionosphere and the interplanetary medium caused these fluctuations. 
Ananthakrishnan & Dennison (1989) saw this effect in an observation of 3C 286. These 
phase fluctuations contain information on the strength and spectral index of electron 
density fluctuations in the solar wind (Cronyn 1972), and will be analyzed in a subsequent 
paper. We removed this effect by subtracting a phase <pbb{t) averaged in frequency over the 
14 MHz bandwidth of our experiment. Since Av < 1 MHz for our sources we expect this 
not to affect the visibility phase fluctuations associated with interstellar scattering. 

The visibility of the pulsar 1919+21 on the 7812 km Owens Valley - Westerbork 
baseline is plotted in Figure 1. The visibility amplitude shows fluctuations in both time 
and frequency characteristic of scattering by the ISM. The visibility phase shows no such 
variations, indicating that the scattering disk of the pulsar is unresolved on this baseline. 
We computed the autocorrelation function of the real and imaginary components of the 
visibility, and used the single-point offset to measure the phase structure function (Gwinn 
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Fig. 2. — The observed FWHM angular diameter 6jj of the average scattering disk plotted 
vs. the FWHM angular diameter 9 U expected for scattering through a uniform medium. If 
a uniform medium is responsible for the scattering of these pulsars, points will lie on the 
dotted line 9 H = 9 U . The dashed lines delimit the range of error in 9 U caused by a factor 
of 2 uncertainty in pulsar distances. The angular diameters of three nearby pulsars are 
inconsistent with a uniform distribution of scattering material. 



et al. (1996)). 

D _ ((Im[V(x,x + B;^t)]) 2 ) 
* 4((Re[V(x,x + B;^,t)]) 2 ) 

The single-point offset of the phase autocorrelation function suppresses noise, but preserves 
the contribution from interstellar scattering because these fluctuations correlate over many 
points in time and frequency. We were unable to detect visibility phase fluctuations for any 
of the 5 pulsars. Table 1 lists upper limits on and 9h and lower limits on d s . 

In Figure 2 we plot 9h vs. 9 u for 11 pulsars. All detections are taken from Gwinn 
et al. (1993) and references therein. The upper limits are from Table 1 except for the 
pulsar 0531+21v. This upper limit represents the angular broadening 9 U associated with 
the long term variable component of the Crab pulsar's temporal broadening r. Isaacman & 
Rankin (1977) found both a variable component and a steady component in the temporal 
broadening of this pulsar. The variable component of r changes on timescales of months 
to years, and is thought to be due to scattering within the Crab supernova remnant. 
The steady component of r arises from scattering in the ISM, and the associated angular 
broadening 9 U is signified by 0531+21c. The measured angular broadening of this pulsar is 
consistent with scattering in the ISM (Vandenberg 1976). 
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6. DISCUSSION 

The scattering properties of most distant pulsars are consistent with scattering through 
a uniform medium. The two exceptions are the young pulsars 0531+21v and 0833-45. 
These pulsars are scattered by their enclosing supernova remnants, and the large distance 
to these remnants produces scattering disks with smaller angular diameters. The scattering 
disks of the three nearby pulsars 0329+54, 1133+16, and 1919+21 are too small to be due 
to scattering in a uniform medium. If the scattering material is concentrated in a thin 
screen, the lower limits on d s for the pulsars 0329+54, 1919+21, and 2016+28 preclude 
scattering from the surface of the Local Bubble. With the exception of 1133+16, all pulsars 
in Table 1 have galactic latitudes |6| < 5°. Like low- latitude extragalactic radio sources 
these pulsars might be scattered by HII regions or supernova remnants in the galactic disk. 
However, it seems unlikely that all of them would be scattered by small features close to 
the pulsars. We can use the lower limit on d s to compute a lower limit on the component of 

— * 

the velocity V SCT of such a feature. This lower limit occurs when the screen is moving in the 
same direction as the pulsar. These lower limits appear in Table 1, and indicate that the 
scattering material must be moving in the same direction as the pulsar. While the required 
velocities for the scattering material are not large, they suggest that the scattering material 
may be associated with the pulsar. 

To explain these results, we argue that material close to the pulsar dominates the 
scattering for nearby pulsars, producing small scattering disks. This scattering could 
be caused by strong electron density fluctuations present at the shock front between an 
electron positron wind and the ambient ISM. Relativistic winds from pulsars can produce 
X-ray synchotron nebulae (Cheng & Helfand 1983, Seward & Wang 1988), indicating the 
possible presence of such a shock front. The pulsar 2224+64 is an example of an older high 
velocity pulsar that has long since left its supernova remnant and whose wind has evacuated 
a cavity in the ISM (Cordes 1993). These relativistic winds may be present at low levels 
in all pulsars, but may nevertheless induce the dominant electron density fluctuations for 
these nearby pulsars. 

For more distant pulsars the integrated interstellar electron density fluctuations 
grow to dominate the scattering, and produce scattering disks consistent with a uniform 
distribution of scattering material. In this model the Crab pulsar would be an intermediate 
case, in which the ISM dominates the angular broadening but the effects of scattering in 
the supernova remnant are apparent in the variable component of the temporal broadening. 

We thank Omer Blaes and Jon Arons for several useful discussions. This work was 
supported in part by the National Science Foundation grant AST 92-17784. 
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Table 1. Scattering Properties of Nearby Pulsars 
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0329+54 


1430 


0.028 


5.7 


4.3 


5730 


<0.21 


<2.6 


>820 


146 


145 


>4 


1133+16 


270 


0.43 


0.37 


2.7 


7870 


<0.17 


<1.5 


>160 


77 


475 


>50 


1919+21 


660 


0.17 


0.92 


2.5 


7810 


<0.085 


<0.77 


>560 


150 


122 


>9 


1929+10 


170 


0.99 


0.16 


2.1 


7940 


<0.30 


<2.7 


>39 


351 


86 




2016+28 


1100 


0.046 


3.5 


3.8 


6510 


<0.28 


<3.1 


>470 


463 


12 





a Distances are from Taylor, Manchester & Lyne (1993). 

b Scintillation bandwidths are from Cordes (1986) and have been scaled to 326 MHz using the 
relation Au oc v AA . 

c Broadening times are calculated using the relation r = 1/2-7tAza 

d Scintillation timescales and proper motion velocities are taken from Gupta (1995). Scintillation 
timescales have been scaled to 326 MHz using the relation t- lss oc u 12 . 
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